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I. INTRODUCTION 


The ccmplexity involved in designing a control system 
cf reasonably low crder for dynamic systems cf substantially 
higher order is a subject about which much has been written. 
The list cf references at the end of this study include only 
a small fertion of the material available concerning this 
subject. A new development has been intrcduced by Maurice 
F. Hutton and Bernard Friedland, [Ref 1], which provides a 
systematic and analytical approach to obtaining 
reduced-crder approximations from higher-crder transfer 
furcticns. This .method is referred. to as the. "Routh 
Approximaticr Methcd," and is based on an expansion that 
uses the Rcuth table of the original higher-order transfer 
function. Tke desire for low ocrder models to simulate 
higher-order practical systems, such as electrical fewer 
plants, chemical processes, aircraft model designs and 
Slectrenic. Circuitry is .obviews when, dealing .with. such 
ccmuplex systezs. The "Routh Approximation Method", is thus, 
an extremely powerful tool to utilize in dealing with these 


types of preoklems. 


In ccmpariscn to other available methods presently 
used, the Routh methcd has many advantages. The widely used 
"Dominant Pcle", approximation method , which is based on 
approximating a system by Wtilezing the poles, or 
characteristic rocts nearest the imaginary axis, has the 
disadvantage that the roots of the characteristic equation 
must be found. For a very~high-order system, this is not a 
trivial task. The "Pad'e" approximation is Fased on setting 
the numerator and denominator crders to a desired value, the 
ccefficierts are then chosen so that the Taylor series 


expansions cf the approximant and the original transfer 
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function agreé€é in aS many terms as possible. This method 
preduces accurate results, however, it is limited in 
applicaticn tc single-input, single-output system analysis. 
In additicn, an unstable approximant may be obtained froma 
stable system since the approximant's poles depend on toth 


tke original eguation's numerator and denominator. 


The "Rcuth Approximation Method," preserves stability 
a 6othe” Criganal transfer’ function’ is stable. It provides an 
efficient means of obtaining lower-order approximants for 
multiple inputs or ottputs and is very adaptable to computer 


programming. 


This thesis is concerned with the development of a 
computer fregram which utilizes the "Routh Approxization 
Method," to obtain lower-crder transfer functions from 
higher-crder functions and compares the resulting equations 
for ene "Cligqiutal eguations “by - graphically displaying the 
response cf each tc various input excitations. The program 
was designed to provide output data which is useful in 
determining which degree of approximant is best suited to 
Simulate tke original higher-order equation. The crders 
available are the first through the fourth order reduced 
tzanster ~ functions. The program is capable of reducing 
transfer functions up to, and including the tenth order and, 
without loss of generality, may be extended to handle an 


almost unlimited order. 


For iliustrative purposes, a seventh crder system was 
used, as an exampie, to demonstrate the computer program's 
capabilities. This particular system was chosen merely to 
demcnstrate the simplicity involved in utilizing the program 


and has ne physical relationship to any real system. 
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Pier sNAtURE CF THE PROBLEM 


In eppreximating higher-crder systems by lower-order 
models, a linear system is desirable since the linear 
characteristic equations are less complicated than the 
non-linear e¢équaticns. Therefore) this study is primarily 
ccncerned with linear, time invariant characteristic 
eguations, since the objective in developing the computer 
prcgram was tc illustrate a simple analytical approach to 
Obtaining lcwer crder characteristic eguations from higher 


order systems. 


Bate is obvicus that a reduced-order model cannot 
characterize a given system as accurately as a higher-order 
model. The validity of the lower-order model is based upon 
its degree of success in approximating the higher order 
systen in representing the characteristics of primary 
interest. 


Interpreting the soluticn of a higher-order systen 
Geren cesults in computational difficulties which are 
recuced ry appropriate selection of a reduced-order 


approximaticn. 


Ideally, a reduced-order model would approximate the 
higher-orcer system in both lcw and high frequency ranges. 
im doing Sc, scme accuracy is lost in compensating fcr the 
different responses of the systen to variations in 


freguency. 
The lcw frequency model more closely approximates the 


higher-order system than a model composed cf both lew and 


high crder frequency characteristics. The procedure 
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utilized in this study places emphasis on the low freguency 


model. 


For unstable systems, the program described herin is 
still valid, however, the original higher-order system must 
first be mcdified by shifting the imaginary axis prior to 
the apprcximation. Thus £0r an unstable transfer 
function,#(s), the equation must be changed tc H(Sta), where 
a is chcser sufficiently large so that H (s) is 
asymptotically stable. This procedure is described in 
decals ft [Ref 1], ,.pp.332. 
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Peis. COMPUTER PROGRAM CRITERIA 


A. GENERAL 


The program, referred to aS ROUTH1, was written with 


the. following criteria: 


ae Minimum utilization of computer time 


ROUTH! consists of less than 150K of storage and takes 
less than 2 minutes of computer time. Phas. net “ony 
conserves efficiency, but also provides the user with the 


desired data at a minimum cost. 


Input data reguired consists of eleven data cards for 
Maximum utilizaticn of the program. Emphasis is placed on 
the ease cf using the program to obtain the desired results 


With minimum time expended on computer programing. 


a. Mew ttrst,, second, third and fourth order 
gpecoximants to the original equation are printed in 
toeansrer: function fornat. Both numerator and denominator 


coefficients are printed in ascending powers of S. 
ihe The roots of the original and reduced equations 


are provided to enable the user to study the response cf the 


systems in the frequency domain. 
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ef three fossible runs. 


di. Ghovce OF VartabLes For graphical output with up 
to four curves may be plotted separately cr ail cn one 


graph. 


e&. Gtaphical output response to input excitations 
consisting Gia Step, kale, Or Sinuscidal input are 
available. The graphs display the original output response 
compared to the lower order response and the error is 
pictted tc display the differences in response. Muitiple 


inputs may be used. 


Be. SUBRCUTINES 


Two subroutines are utilized in ROUTH1. The rocts of 
the original and reduced equations are determined by 
Subroutine PROD , which was taken from [Ref 3], and the 
tables and graphs are determined by subroutine REDUCT1 which 
fad) weditication Gf) INTEG? ,»1Ref 3]. 
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IV- FOURTH ORDER EXAMPLE 


A. THREE STEP PROCEDURE 


To, aidwustratce, the “Routh. Approximation Method”, A 


fourth order example was chosen for simplicity. 


The method, Gescribed in detail in [Ref 1], consists of 


three basic steps. The following transfer function is used 
to illustrate the procedure: 


H(s) = Lt 


20 + 328 + 248+ 89+ ¢ 


The first step is to compute what are termed Alpha and 


Beta coefficients from the Routh Table shown on page 19. 


20 24 He 


Alpha, =0.625 19 1 
Alpha,=1.684 aes 
Alpha,=3.008 1 

Alpha ,=6.3158 


Beta 


Beta, =0.03125 =O) AS 
Beta,=0.0 07.0 
Beta,=0.03958 
Beta,=0.0 


16 


osHUOF 


TIUCEIORS Were wae 


sig stensegiliy’ Gig 
(i> €@sv slqmege tebao me 


- sje. 6 wditogeb ,badtvonr ef] 


13 ‘orth st eff Deeete sigsdo® 
~ +: sxebarergegy aiid StaitsaL ¥ 


— \e {818 


~ scugnes of ek gage Sestiiam 
q ye at; en cust esnsistt rete 


r; 


Mm 


=- 


Ss hs 
f ad 


ee SE A ae ae 
. 
a 


Step two in the procedure is to obtain what are termed 
the Routh convergents, which are based on the following: 
Letting A, (s) and B (s) denote the denominator and the 


numerator, respectively, of the eon Routh convergent, i.e. , 


A, (Ss) = Alpha,S + 1 

B, (s) = Beta, 

A,(s) = Alpha,Alpha,s* + Alpha,s + 1 
B. (s) = Alpha,Beta,S ~ Beta, 


The general expression from [Ref 1] is the following: 


A, (Ss) = Alpha, (s)A,_, (s) + A.» (s) 
By. (s) = Alpha, SB, _, (s) + By._> (Ss) + Beta, KU eyo. <e 
with A_, (s)=0.0 B_, (s)=0.0 
Ay (s) =1.0 By (s) =0.0 


the Routh CoOnvergencts for the fourth order example are 


the following: 
R, (s)= G-03125 
O62 552+ a 


R,(s)=__ 0.056316 
1.05263S° + 1.68425 + 1.0 
2 


R,(s)=__0.1588335" — 0.0083 
Beicgse = Sedcec- § 3.6335 + 1 
a 3 


20s? + 328° + 248° + 8S + 1.0 
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The third, and final step in the procedure is to apply 


what is termed a reciprocal transformation, defined by 


ak 
H. (s) = sxR, (1/s) 
which is merely a reversal of the order of the polynomial 
coefficients. Thus, for the example given, the reduced order 


approximations are given by the following: 


H, (s) = 0.03125 


S + 0.625 
H,(s) = 0.056316 
S-pieeedos @ Vs05063 
H3(s) = -0.0083S* + 0.158333 
s? + 3.6338 + 5.0665 + 3.166 
ots) = “i 
4 ri 3 5 


Spot oom tess + eo2s + 20 
As expected, the fourth-order approximation is the same 
as the original equation. 
The poles of the approximants are illustrated in table 
IV.2. The poles of the approximants approach the dominant 
poles of the original higher-order equation as the order of 


the approximant is increased. 
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Table IV.2 
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Ve. COMPUTER EXAMPLE 


he S@eVENTH ORDER SYSTEM 


A seventh order system with the transfer function given 


by [Ret 2:], as 


384x10/ 


3 2 


Be 7 % z r 
s/+4325°+62670S°+3615900$~*+75114000$2+553920000S 
7 


- -+1443200000S+384x10 


or in factored form as 


384x10/ 


G(s)= 
($*+2S+10) (S+10) (S+20) (S+80) (S+120) (S+200) 


was reduced to the low order approximants by 
utilizing ROUTH1 and the response to input excitations, 
consisting of a Step and Ramp are illustrated in figures 5.1 


througa 5.6. 


The rocts of the system and its low order appreximants 


are illustrated in table V.1. 
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Figure 5.6 


Vio COMPARTSON TO OTHER METHCDS 


A. DOMINANT FOLE METHOD 


The Dominant Pole Approximation method is based upon 
Gerla zZing the poles closest to the imaginary axis. fhe 
equation must. be factored to obtain the characteristic 
rocts. Fcr the given seventh-order system, the lower- order 


eguations are given by this method as the following: 


10 


2 


H(s) =——————_—— 
S +2S+10 


100 


3 Z 


H(s) =—————_>_- 
S~+14S°+30S+100 


2000 


4 S 2 


H(s) =-—————— 
S 845) 3210S) +700S+2000 


Graphical flots of the Dominant Pole reduced 
SqUuatitOns, in fesponse to Step inputs, are illustrated in 
maguees 6.1 through 6.3 in comparison to the original 
equation and figures 6.4 through 6.6 illustrate the 
ccmpariscn to the Routh equations. The error between the 
systemsomas also plotted. Table VI.1 gives the analytical 


data for perfcrmance measure comparisons. 
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Be SLeERATEVE OPTEALZATION METHOD 


The fcllowing equations, representing the seventh-order 
were taken from [Ref 2]. This method makes use of 


example, 
an iterative tinimization technigue to locate the best pole 
and zero locations for the lower-order nodels. 
7.203856 
LC eae ae 
So +l. Jseel6s+7. 2033856 
52.2901 
li: \S |= Se eae ya 
Saae LoOoSs tL 92 50LS5S+5 25286) 
1470.1403 
Bg) eee 
on eile O20 4S504+209.3842S +552.55241S+1470 21403 
Plots of these equations versus the Original 


and the Routh equations are illustrated in figures 


eguation 
imp tS Table Viet 


6.7 through 6.12 in response to Step 
provides numerical data for performance measure comparisons. 


30 


a : % 


COW?GN WOTTSETETEGO SVEPARSES” 


; ra ion ie 

793 ,2celPapee pniwchios sag J 
¢ Jes ce gece? s76u . apde 
, eae 


pisdoe? aolsémteatate sviseaem 


' . « 
wool edi IoF sfo.ssoee Ose 


sPacos. 7 


"> ku, ee —— a 


2 Sees PSP OL Oke ReOke, BSS a | 
Be 

‘fe 

Jsuye @60a2 Fo asoit 


5 


ecic C22 SAEOLAERRS stot was bas 6m 
nt gs4t oF sunoqewa) GES? ye epee 

ni - 

sos soreesodaeq 569 tab Daeweaem 


x 
a 


alee 


Gy PRESENTATION OF DATA 


Comparisons were made between the "Routh Approximation 
Method", and that of the Dominant Pole and Iterative 


Optimization methods previously described. 


Therphasis £On) COMbperison consists, ©@f the following 


eriteria: 3 


1. Peak Oe ee ax 


ss 

ioe OTS ILERIA Rabe inley— T,=time fOr -Y (E)--to-reach-O~5 vee 
the very first time. 

S14 RGLSIG) Genus se T =time Ore) ete. gorccom O.i “tS 
O29 Gf the final value. TW=1/8W. 

a. | Settling Time--~-T_=time at which Ge YS 


5. Graphical Representation in response to Step inputs. 


6. Average Error----J= JZE/t,| ; t,=integration steps 


Table VI.1 illustrates the above comparisons. 
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SYSTEM | 
TYPE Mpt Td me Ts g 


" 
o-4 | 1.325 |o.s33 lo.7a 5.866 \sx10 : 
i i 
i 
D-4 115338 10.528 
| JS eee 
R-3 e247 lie 533 
| | : 
O-3 | #. 320 O25 33 ‘ 
! 
D=3 | 1.346 [0.502 
i pee ae eee es Lee 
R-2 | 1.083 10.400 jo.712 5.966 0.062 
1 } | : he one ew 
i { | 
027) ey i lle CRE ao7it 5.866 lo.037 : 
= Oo ea | 
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Table 


Original Seventh Order 
Symbols Fd 


Routh Aporoximation 
Dominant Pole Approximation 
Optimum Minimization Method 
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VII. CONCLUSION 


The “graphical data and tabulated analysis of this 
thesis irdicate that the Routh Approximaticn Methecd is a 
valuatle fermulation technique which produces very 
satisfactcry results in acquiring approximations to 


higher-order systems with minimum cost. 


In ccmpériscn to the other methods cf analysis, the 
Rcuth Appreximation Method offers a quick and easy 
analytical approach to obtaining lcw order models. Unlike 
the Fad'te Approximation, the Routh method ensures statility 
of the lcwer-crder models, if the original higher-crder 
system is statle. The original system need not be factored, 


as in the Dcninant Pcle Approximation method. 


They GCcaputer program, ROUTH1, utilized to acquire the 
reduced crder equations, the roots cf the eguations and the 
Gaaphical. pricts. and numerical tables, takes considerably 
Ressa tame, than the minimization technique utilized in 
reference 2. However, the minimization technique operates 
efficiently without prior knowledge of the system's transfer 
function. in. the. Routh. -metaodthe...availability.c£ the 


transfer function is a necessary requirement. 


In ccmparing the low order equations of the various 
methcds discussed, the Routh method has proven to be a valid 
and efficient solution to the problem of cbtaining gocd low 


crder approxigants to complex higher-order systems. 
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